Impact Parameter Dependence in the 
Balitsky-Kovchegov Equation 



Takashi Ikeda a , Larry McLerran a ' b 

a RIKEN BNL Research Center, Brookhaven National Laboratory, Upton, NY 

11973 USA 

^Physics Department, Brookhaven National Laboratory, Upton, NY 11973 USA 



Abstract 

We study the impact parameter dependence of solutions to the Balitsky-Kovchegov 
(BK) equation. We argue that if the kernel of the BK integral equation is regulated 
to cutoff infrared singularities, then it can be approximated by an equation without 
diffusion in impact parameter. For some purposes, when momentum scales large 
compared to Aqcd are probed, the kernel may be approximated as massless. In 
particular, we find that the Froissart bound limit is saturated for physical initial 
conditions and seem to be independent of the cutoff so long as the cutoff is suffi- 
ciently large compared to the momentum scale associated with the large distance 
falloff of the impact parameter distribution. 



1 Introduction 



The issue of impact parameter dependence of parton distribution functions is 
an old one. Recently, it has been possible to address this issue quantitatively 
in the context of the BFKL and the Balitsky-Kovchegov evolution equations 
[l]-[3]. One would like to compute the parton phase space density, 

dN gluon 1 2N C 0(y,fc T ,6) m 
dYd 2 k T d 2 b a 5 (27r) 2 7r 2 k\ 1 ' 

where Y is the gluon rapidity, kj- is its transverse momentum and b is the 
transverse spatial point where the gluon density is measured. The quantity <j) 
is related to its coordinate space analog by 

N Y (r T , b) r d 2 k T jk t .? t <t>(y,k T ,b) 

4 J (27T)2 e k\ { ' 
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In the Color Glass Condensate description of high density gluonic matter, 
valid at small x [4]- [7], 



N Y (f T ,b) = ^< Tr{\ - U\x)U(y)) > (3) 



where U is a line integral of the gluon field which goes over the rapidities of 
the gluons, and where 

f T = x T - y T (4) 

and 

g + ( 5 ) 
2 v 1 

As r T — > 0, TV — > 0, and as ty — > oo, TV — > 1. This follows because at large 

distances, the expectation values of < U(x)U(y) >^< U(x) >< U(y) >, and 

< U(x) >^ by gauge invariance. This means that correlations should fall 

off at long distances, and it means that isolated source of color charge cannot 

propagate in QCD. 

The distribution function N has been argued to satisfy a non-linear evolution 
equation, the Balitsky-Kovchegov (BK) equation. It is usually written in terms 
of the variables x, y rather than a relative coordinate and impact parameter 
r, b. (Below, we drop the arrow and subscript T representing the transverse 
variables.) 

AiW (x, y) nj^ u -IL _ - (N Y (x, Z ) + N Y {z, y) 

-N Y (x,y)-N Y (x,z)N Y (z,y)) (6) 

where 

a = ^ (7) 

Let us see what happens if we try to find a solution to this equation which 
has exponentially falling behavior at large impact parameter. Let us take as 
initial condition for this equation [8] 

N (r,b) = l-exp(-f(r)e- 2 » b ) (8) 

In this equation, ji should be taken to be m n for the long distance falloff, since 
this should be controlled by isosinglet exchange, and the 2 pion state is the 
lowest energy strongly interacting state with these quantum numbers. 

At large impact parameters, for fixed r, this goes into the factorized form 

hmiV (r,6) = /(r) e - 2 ^ (9) 
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On the other hand, if f(r) increases as r increases, then for some large enough 
r, the exponential on the right hand side of Eq. (8) becomes small, and iV — > 
1, its saturated value. This is the behaviour one expects from the operator 
definition of the distribution function N, since we expect that 



should approach zero at large r. This requires that /(r) be singular at large r. 
It is also true that if we make / grow too rapidly as r — > oo, we may introduce 
a stronger dependence on the infrared than desired. In what follows, we will 
choose initial conditions so that / goes not faster than r 2 as r — > oo. This 
constraint on the growth at large r is not spoiled by evolution, as we will see 
numerically below. Specifically, we shall choose 



Here, c is the mass- dimensional constant corresponding to the initial satura- 
tion scale at zero impact parameter. 

Let us see if the generic features of this form can survive evolution in rapidity, 
according to the BK equation. Let us consider Ny{x,y) for x = R and y = 
where R ^> l//i. After one iteration in time (rapidity), the contribution 
N Y {z,y) in the BK integral is big for z ~ 0. The kernel is of order I. This 
generates a non-exponentially falling contribution. Next if we iterate again, 
and look at Ny{x,y) for x~y~R(r~l/[/, and b ~ 2R), we see that 
iV ~ 1/R A . This means that the solution evolves towards a power law falloff 
in the impact parameter, not an exponential as one would have thought from 
the initial conditions. 

The basic reason that the BK equation generates a bad behaviour in impact 
parameter is because of the massless nature of the kernel of the BK equation. 
The point of this paper will be to show that if one regulates the kernel by 
introducing an infrared cutoff, then the exponential falling nature of BFKL 
which may be built into an initial condition will be maintained by evolution. 

What we will argue is that the BK equation may be replaced by an equation 
which has no diffusion in impact parameter. This equation is well defined with 
no infrared cutoff in the kernel. The corrections to this equation which involve 
diffusion remain small so long as the infrared cutoff is larger than or equal to 
fi/2, the scale associated with the impact parameter falloff and so long as one 
measures the distribution functions at momentum scales Q 2 ^> A.q CD at not 
too large impact parameters. Roughly speaking, Q ~ 1/r, so one probes the 
short distance structure of /(r, b). 

This form of the evolution equation is sufficient to establish the Froissart 
bound for high Q 2 processes [9]- [11]. It is reasonable to assume that this bound 



< U(r + R)W(R) > 



(10) 



f(r) 




(11) 



1 + or 
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is also true for total hadronic cross sections where Q 2 ~ Ag CD . We find that 
the Froissart bound is saturated, and the scale associated with the cross section 
reflects the initial conditions, not the scale in the cutoff, so long as the scale in 
the cutoff is greater than or equal to one half that scale of falloff in the initial 
impact parameter profile. 



2 Regularizing the BK Kernel 



Let us try the simplest possible modification of the BK kernel in Eq. (6). We 
define 

Here, the parameter A represents the non-perturbative mass scale. We will take 
A as a parameter in our analysis below. The smallest physically acceptable 
value for A we believe is A = m^. It might be larger, since in the evolution 
of the distribution function, the growth of distribution functions might be 
associated with vector meson exchange, as is true in some models. The generic 
conclusions we will draw below rely more on A being larger than ~ m n , than 
on its specific value. In the limit that A — > 0, this is the kernel of the BK 
equation in Eq. (6). 

To understand how the infrared cutoff in the kernel above preserves the fea- 
tures of our ansatz for the initial condition for N, we write the BK equation 
with the modified kernel in terms of relative coordinates 



^N Y (r,b) = aJ ^ |i _ z| ;]' i + z|a e-^'e-^l 

-N Y (r, b) -N Y ( r --z, l -{z+ r -)+ b)Ny(^ + z, \{z -£) + &))■ (13) 



Let us check how such an ansatz for the kernel affects the solution. We assume 
the initial condition is given by the ansatz of Eq. (8). 

First consider what happens if we keep r finite and go to large b. Unless 
z ~ ±26, the functions N in the BK equation are exponentially small. If this 
is satisfied, then the kernel of the equation is of order e~ 4Af> so that we pre- 
serve the form of the solution so long as A > /x/2 because the distribution in 
the kernel is factorized in Eq. (9). (We expect that ji > m n for the large dis- 
tance asymptotics of the regulated BFKL kernel, as noted in Sec. 1. Therefore 
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integral 




Fig. 1. Impact parameter dependences of (14) for various A are plotted for r = 0.1 
at cos# = 0. Solid, short dashed, dot-dashed, dashed and dotted lines correspond 
to those for A = 0, 0.1, 0.25, 0.5 and 1. 

this constraint is a rather weak one.) In order to be precise, we evaluate the 
following integral corresponding to the right-hand side of the BK equation 



' ( x V) 2 e - X \x-z\ e -X\z-y\ 

(x — z) 2 (z — y) 2 

x (n (x, z) + N (z, y) - N {x, y) - N (x, z)N {z, y)) , (14) 



with the initial condition (8) together with (11). Throughout this paper, we use 
the unit /i = 1 and set c = 8.0 in (11) in numerical calculations. Figure 1 shows 
the impact parameter dependence of this integral for dipole size r = 0.1 and 
orientation cos 8 = f ■ b/rb = 0. Here, we have assumed cylindrical symmetry 
and dropped the dependence on the azimuthal angle. This clearly shows that, 
while the power law tail is generated for A = 0, the exponential falling tail 
in the initial condition can be maintained for A > fi/2. This power law tail 
can be seen in Fig. 7 of [15] where the same integral with no infrared cutoff 
has been evaluated for the initial condition of the Glauber-Muller form with a 
steeply falling profile in b, that is, e~ b2 . Therefore, the massless property of the 
BK kernel, not the form of the (initial) distribution, may originate the power 
law tail in impact parameter. For A < /i/2, the behaviour of the integral can 
be shown to be e~ 4A6 /6 4 , from the structure of the integral, and this fits the 
observed distribution well. 

Now consider what happens if we take the limit that r — > oo with b ^> l//i, 
but b fixed. In this case, the kernel vanishes, and the distribution function 
maintains its original form, as seen in numerical solutions below. This is good 
because in this region there is little matter, but one is in the non-perturbative 
region. On general grounds, the distribution function should tend to one in 
this region, and should not be much affected by evolution due to the presence 
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of the Color Glass Condensate. 

Note also that if we look at the problem term which generated long distance 
singular behaviour in our original analysis of the massless unregulated kernel, 
we need to look at r ~ R and b ~ R/2. The potentially dangerous region 
comes from z ~ —R/2. Then this contribution is suppressed by e~ XR from the 
regulated kernel which is small so long as A > /i/2 as was the case above. 

To summarize, the modification of the kernel we propose does in fact have 
the correct properties to preserve the behavior we expect of N Y (r, b) based on 
general principles. 



3 An Equation Ignoring Impact Parameter Diffusion 

3. 1 Discussion 

The Balitsky-Kovchegov equation as written in the r, b basis as in Eq. (13) 
generates a power law tail in inverse impact parameter if the kernel is not 
regulated. On the other hand, the regulated kernel generates no such tail. 
This suggests that if we compute the distribution function N(r, b) on scales 
of r which are small compared to the scale of variation in b, then one may be 
able to ignore impact parameter diffusion. 

Looking at Eq. (13), we see that if the exponential tail of the impact parameter 
cutoff satisfies A > /i/2, then the dominant contribution in the integration over 
z comes from z ~ r. In the region where z ~ l//i, we get a contribution of 
order r 2 /i 2 . Clearly, ignoring impact parameter diffusion should be valid so 
long as we have a cutoff kernel. 

In practical terms, this means that a good lowest order approximation should 
be generated if we expand the terms inside Eq. (13) 

N(r',b')=N(r',b + Ab) (15) 

in powers of Ab, and keep the first non-leading terms. The lowest order term 
in this gives an equation local in impact parameter (after transforming z — > 
z-±r)[12] 

— N Y (r,b)=a J d 2 z^-—^-^e- x ^e- x ^(N Y (r - z,b) 

+N Y (z, b) - N Y (r, b) — Ny(r — z, b)N Y (z, b)) . (16) 
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This equation clearly preserves the initial conditions at large impact parame- 
ter. If r <C 1///, it also has the property that one can ignore the A dependence 
in the evolution equation, so long as A > /i/2, as we shall assume. This means 
the solution to this equation has all the properties assumed in Ref. [11], where 
the Froissart bound was argued to be true when measurements were made at 
scales Q 2 » fj 2 . 

The solution to this equation should be roughly of the following form: In 
the limit where the distribution function is itself small and where r <C 1/A, 
the solution is the ordinary BFKL solution of the linear BFKL evolution 
equation times an impact parameter profile. This solution will hold until the 
function N becomes of order 1. This should occur roughly when 

(rQ sat (Y)ye-» b ~l, (17) 



or 

r ~oZm (18) 

where Q sa t(X) is the saturation momentum of the ordinary BK equation. 
For rapidities beyond the critical rapidity for which the amplitude Ny(r,b) 
first becomes of order 1, the amplitude remains of order 1 for larger rapidity 
because the BK equation doesn't violate the unitarity. 

On the other hand, when r > 1/A, the kernel of the BK equation introduced 
above cuts off further evolution, and the function N Y (r, b) is frozen to its initial 
condition. 

It is difficult to get a direct computation of the correction to the zero impact 
parameter diffusion approximation. Naive expansion in a Taylors series around 
zero impact parameter diffusion leads to an equation with mild infrared singu- 
larities. Nevertheless, the estimate made above, that these corrections should 
be of order r 2 /i 2 , up to logarithms, is still true. This of course requires that 
A > /i/2. If this was not true, the tail of the distribution would be modified 
in an unmanageable way. Nevertheless, the leading order solution is valid so 
long as we study small r, and does not requires that A > /i/2. On the other 
hand, at large r, our form of the kernel has the physically plausible behavior 
that iVy(r, b) does not evolve in Y. 



3.2 Numerical Solutions 



In this section, numerical solutions of Eq. (16) for A = and 1 are calculated. 
We used the same parameters as = 0.2, N c = 3 and c = 8 in the initial 
condition (8) together with (11), as in Sec. 2. Since Eq. (16) is local in b and 
our initial condition depends only on r(= |f|) and b(= \b\) but not on the angle 
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Fig. 2. The r dependence of distribution functions for b = 0.13 (left panel) and 4.95 
(right panel). Dashed lines represent the initial distribution. Solid and dotted lines 
correspond to the distributions for A = and 1 at Y = 40. 

8 between r and b, the solution Ny at any rapidity has the same property, 
that is, N Y (r,b) = N Y (r,b). 

The numerical method to solve Eq. (16) is similar to the one for the solutions 
of the full BK equation [15], and highest rapidity in our calculations is 40. 



3.2.1 The dipole size dependence of the solution 

First, let us examine the rapidity evolution of the distribution function for 
small dipole sizes r < 1. In Fig. 2, the distribution functions for A = and 
1 are plotted as a function of r. Left and right panels correspond to the 
distributions at b = 0.13 and 4.95. In each panel, the dashed line is the initial 
distribution, and solid and dotted lines are the distributions for A = and 1 at 
highest rapidity Y = 40. While the distribution functions for A = and 1 are 
in good agreement for small b(= 0.13) as shown in the left panel, the evolution 
for A = 1 is slightly suppressed for large b(= 4.95) in the right panel. We also 
found the exponent for small r changes from r 2 in the initial distribution to 
r 270 with 7q ~ 0.6 — 0.7 at high rapidities, independently of b and A [13]-[14]. 
This 70 is called anomalous dimension and the value is consistent with other 
studies [11]. 

Next, let us examine how the regulated kernel affects the evolution of large 
dipole size. Figure 3 shows the distribution function at b = 4.95 for A = 
1 including the region r > 1/A where the regulated kernel may affect the 
evolution. It is observed in this figure that the evolution of the distribution 
function is suppressed around r ~ 1(= 1/A) and the distribution is frozen to 
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Fig. 3. The distribution function at b = 4.95 for A = 1 is plotted as a function of r 
at various rapidities. 



its initial value in the region r > 1, as we mentioned in the previous section. 
If the initial distribution reaches to its saturated value in the region r > 1/A, 
it doesn't evolve any further in such region and the regulated kernel doesn't 
affect its evolution. This is the case for small b in our initial condition, as you 
can see in the left panel of Fig. 2. 

That said, the dip shown in Fig. 3 for large r is in an unphysical region where 
the separation of the sources of color charge within the dipole are at unphys- 
ical size scales. This is a region dominated by non-perturbative physics, and 
ascribing physics to this region is no doubt problematic. This region causes no 
problem in the solution of the equation we consider, and the dip is presum- 
ably an artifact of an incomplete non-perturbative treatment of the correlation 
function. 



3.2.2 The impact parameter dependence of the solution 



In Fig. 4, the impact parameter dependences of the distributions at fixed r 
are plotted for Y = and 40. It is clearly observed that the exponential 
falling behavior in the initial condition is maintained at high rapidities both 
for A = and 1, as we expected above. While, for small dipole size r = 10~ 10 , 
the distribution functions for A = and 1 at high rapidities are in good 
agreement, the evolution for A = 1 is a little suppressed compared to that for 
A = for r = 0.1 due to the regulated kernel. 
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Fig. 4. The distribution functions for fixed r at Y = and 40 are plotted as a 
function of b. Left and right panels correspond to the distributions for r = 10 -10 and 
0.1. Dashed lines represent the initial distribution. Solid and dotted lines correspond 
to the distributions for A = and 1. 
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Fig. 5. The saturation scale for A = (solid line) and 1 (dotted line) are plotted as 
a function of rapidity. Upper and lower lines correspond to Q S at(Y,b) for b = 0.13 
and 4.95. 

3.2.3 Saturation scale 

Next, let us examine the rapidity evolution of the saturation scale Q sa t(Y)- 
We define Q sa t(Y, b) as an inverse of the lowest r satisfying 

N Y (r = l/Q aat (Y,b),b)=K, (19) 

where k is order of unity. The explicit value of k doesn't matter and k is set 
to 1/2 below. This definition is the same as in Ref. [15]. 

Figure 5 shows the saturation scale for b = 0.13 and 4.95 as a function of 
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rapidity. As you can see, the rapidity evolutions with A = and 1 almost 
coincide for small b (=0.13 in this figure), in other words, the infrared cut- 
off doesn't change the evolution of the saturation scale for the small impact 
parameter. 

For large b(= 4.95), the saturation scale for A = evolves similarly to that for 
b = 0.13. On the other hand, the Q sa t for A = 1 stays in its initial value until 
the rapidity goes beyond a certain value, and then starts evolving. Looking 
carefully at Fig. 3, we can understand this behavior. As rapidity increases, 
the height of the bump of the distribution for A = 1 around r ~ 1 also 
increases. Until the height of the bump reaches to k in (19), the Q sa t(Y, b) is 
frozen to its initial value. (We should note that nothing particularly exciting 
is happening with the correlation function itself when Q sat makes a jump, and 
the jump is no doubt an artifact of its definition, rather than corresponding 
to some rapid change in a physical quantity.) When the height reaches to k 
(at Y ~ 28 for b = 4.95), the Q sat jumps to the lowest r satisfying (19) from 
the initial value l/Q sa t(Y = 0) and then evolves normally because this is the 
region Q sa t(Y, b) > A, e.g., r < 1/A and the evolution of the distribution is not 
affected by the regulated kernel. The point at which the Q sa t(Y, b) jumps from 
the initial value and starts evolving is a monotonically increasing function of 
b. For example, the Q sat at b = 3.3(4.95,7.4) has a jump at Y ~ 20(28,40), 
and for b > 7.4, the Q sat for A = 1 is still frozen to its initial value because of 
the highest rapidity Y = 40 in our calculations. 

It is noted that the slopes at high rapidities are almost same independently of 
b and A, which means the saturation scale at high rapidities can be factorized 
in the form 

Q sat {Y, b) ~ exp{2n\ s aY) x N(b) , (20) 

where N(b) is an impact parameter profile function. In our calculation, X s is 
estimated around 2.1 which is consistent with that from the full BK equation 
[15]. This might mean that the rapidity evolution of Q sa t{Y, b) is not sensitive 
to the impact parameter diffusion. 

In Fig. 6, the b dependence of the saturation scale for A = and 1 is plotted at 
Y = and 40. For small b, especially b < 1, the impact parameter dependences 
for A = and 1 agree well. This is because the evolution of the distribution 
function for small b doesn't change by the vacuum property of the regulated 
kernel, as seen in the left panel of Fig. 2. 

We also found the exponential tail of the saturation scale for A = is main- 
tained at high rapidities, e.g., e~ 2b . Of course, this exponential tail is totally 
different from that from the full BK equation without infrared cutoff of the 
kernel, which reads the power-like tail Q sa t{Y, b) ~ 1/b 1 with 7 = 1.6 — 2.0 at 
high rapidities [15]. The exponential tail is preserved only when the kernel of 
the BFKL equation is properly regulated in the infrared. 
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Fig. 6. The impact parameter dependence of the saturation scale for A = (solid 
line) and 1 (dotted line) is plotted at Y = and 40. 

4 The Growth of Black disc radius and The Froissart bound 



One of the most important issues in the saturation physics is how fast the black 
disc radius in which Ny saturates grows in rapidity. This is closely related to 
the problem of the Froissart bound [9] which tells us that the total hadronic 
cross section is bounded by 

a = ^\n\s/s ), (21) 
ml 

with energy s ~ e Y . We define the black disc radius Rbd^, Y) as 

N Y (r,b = R BD {r,Y)) = K, (22) 

where k is order of unity and k — 1/2 is chosen as in the case of the saturation 
scale. 



In the saturation regime, the total dipole-nucleus cross section is given by 
integrating the scattering probability Ny over the impact parameter, 

a(r, Y) = 2j d 2 bN Y (r, b) ~ 2irR 2 BD (r, Y) . (23) 

If this Rbd grows at most linearly with the rapidity Y, the Froissart bound is 
saturated. Let us discuss this below. 

Figure 7 shows the rapidity dependence of the black disc radius for various 
dipole sizes. As we can see, in both cases for A = 1 and 0, Rbd grows linearly 
as a function of rapidity. The slope at high rapidities for A = is slightly 
larger than that for A = 1 which doesn't depend on the dipole size if it is 
small enough (r < 1). One might think it's a natural consequence of the local 
approximation ignoring the impact parameter diffusion. The important point, 
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Fig. 7. The black disc radius for A = (solid lines) and 1 (dotted lines) is plotted 
as a function of rapidity at r = 0.1,0.01 and 0.001 from the left to right. 

however, is that this linear growth of the black disc radius in Y is expected in 
the full BK evolution equation if the kernel of the equation is regulated by the 
infrared cutoff, although the naive BK equation shows the exponential growth 
and violates the Froissart bound. 



5 Summary 

The origin of the infrared cutoff in our analysis needs some understanding 
beyond what we present here. It is clear that the lower bound we have is a 
weak one and that the generic feature of Froissart bound saturation appears 
to be insensitive to details of this cutoff. It is however not clear to us what 
sets the scale of this cutoff. Is it 2m n , or 2m p ? Might one derive an effective 
meson theory which can generate a BFKL like equation to describe this non- 
perturbative region? Of course, it is purely non-perturbative and cannot be 
explained by the perturbative QCD. 

Although the considerations presented here are far from rigorous, it seems that 
one does in fact have a simple intuitive picture of how the Froissart bound 
arises. It is simply the trade off between the exponential growth of the gluons 
as a function of rapidity filling up and exponentially falling tail of an impact 
parameter distribution set by initial conditions. 
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